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ABSTRACT 

The Milagrito water Cherenkov detector near Los Alamos, New Mexico, 
operated as a sky monitor at energies of a few TeV between February 1997 and 
May 1998 including the period of the strong, long-lasting 1997 flare of Markarian 
501. Milagrito served as a test run for the full Milagro detector. An event excess 
with a significance of 3.7 a from Markarian 501 was observed, in agreement with 
expectations. 

Subject headings: BL Lacertae objects: individual (Markarian 501) - 
gamma-rays: observations 
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1. Introduction 

Very High Energy (VHE) 7-ray astronomy studies the sky at energies above 100 GeV. 
To date, 4 galactic and 4 extragalactic sources have been identified as VHE sources (see 
Ong 1998 and [Hoffman et al. 1999| for recent reviews). The four extragalactic sources, 



Markarian 421 (Mrk421, z = 0.031) ([Punch et al. 1992|) , MrkSOl (z = 0.034) ( IQuinn et| 



|al. 199q ), lES 2344+514 (z = 0.044) QCatanese et al. 1998| ) and PKS 2155-304 {z = 0.117) 



( IChadwick et al. 1999|) are relatively nearby objects of the BLLac subclass of active 
galactic nuclei. A characteristic feature of BL Lac objects is their rapid flux variability at 
all wavelengths. Flaring activity at TeV energies has been observed both from Mrk421 and 
Mrk501, ranging from variability times scales of minutes ( paidos et al. 199€ ) to months 
( IFrotheroe et al. 1997|) . 

Source detections and analyses at VHE energies are currently dominated by the highly 
successful atmospheric Cherenkov technique. Cherenkov telescopes are excellent tools for 
the detailed study of point sources and their sensitivity has been significantly improved 
over the past few years. Strong fiaring activity of Mrk421 and Mrk501 can be detected 
with less than an hour of observation time per night. 

To complement the pointed atmospheric Cherenkov telescopes, there is a strong case 
for wide-aperture instruments monitoring the sky with a high duty cycle and performing 
an unbiased search for new sources and source classes. The price to pay for overcoming 
the limitations of atmospheric Cherenkov telescopes is a loss in sensitivity for individual 
sources. To date, no unambiguous detection of a steady TeV source has been established 
with an air shower detector. 



The Milagro water Cherenkov detector (|McCullough et al. 1999D near Los Alamos 



New Mexico, at latitude 35.9° N and longitude 106.7° W, a first-generation all-sky monitor 
operating with an effective energy threshold below ITeV, started data taking in 1999. 
Milagrito ( [Atkins et al. 1999] ), a smaller, less sensitive prototype of the top layer of Milagro, 



took data between February 1997 and May 1998. Milagrito was located at the same site and 
served mainly as a test run for studying specific design questions for the Milagro detector. 
Nevertheless, Milagrito operated as a fully functioning detector and took data during the 
strong, long-lasting flare of Mrk501 in 1997. During this flare, Mrk501 was intensively 
studied with several atmospheric Cherenkov telescopes (see [Protheroe et al. 1997[ ). Detailed 



flux and spectral studies have been pubhshed from data taken by the Whipple telescope 
on Mt. Hopkins (Arizona) between February and June 1997 ( ^amuelson et al. 1998| ), and 
the HEGRA stereo system of Cherenkov telescopes on La Palma (Canary Islands) between 
March and October 1997 ( [Aharonian et al. 1999 ). Although they do not cover the same 



observation times, the average fluxes measured by Whipple and HEGRA agree extremely 
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well in both shape and magnitude, and they both indicate an energy spectrum that deviates 
significantly from a simple power law. Using an average fiux as measured by Whipple, 

. s -2.20 ±0.04 ±0.05- (0.45 ±0.07) logio (£^/TeV) 

JJE) = (8.6 ± 0.3 ± 0.7) X 10"^ I — — ) m"^ g-i TeV"\ 

Kiev J 

(1) 

simulations suggest that the observation of a statistically significant excess from MrkSOl is 
within the reach of Milagrito. 

Observations with atmospheric Cherenkov telescopes do not cover the time between 
October 1997 and February 1998, when MrkSOl was visible only during the day time. 
When atmospheric Cherenkov telescopes resumed observations in 1998, they observed a 
relatively high fiux for a few days at the beginning of March, but the fiux quickly decreased 
and was considerably lower for the rest of 1998 than in 1997 (IQuinn et al. 19991) . As an 
instrument that is insensitive to sunlight, Milagrito continued to monitor MrkSOl in late 
1997 and early 1998. In this Letter, we present the results of an analysis of Milagrito data 
on MrkSOl. 



2. The Milagrito Detector 

As source spectra tend to be falling power laws, a large detector area is essential for 
a sufficient rate from point sources in the VHE region. The restriction to earth-bound 
detectors makes the detection of the primary 7-rays considerably more complicated, as 
the primary particle generates a cascade of secondary particles in the atmosphere, an "air 
shower." Air shower detectors have to reconstruct the properties of the primary 7-ray from 
the secondary particles reaching the detector level, and any 7-ray signal has to be observed 
in the presence of a large isotropic background from cosmic rays. To achieve sufficient 
sensitivity at TeV energies, a high altitude location and the ability to detect a large fraction 
of particles falling within the detector area are crucial. 

Milagrito was a water Cherenkov detector of size 3S m x SS m x 2 m, located at 26S0 m 
above sea level (7S0gcm~^ atmospheric overburden) in the Jemez Mountains near Los 
Alamos, New Mexico. The project took advantage of an existing man-made rectangular 21 
million liter pond. A layer of 228 submerged photomultiplier tubes on a 2.8 m x 2.8 m grid 
detected the Cherenkov light produced by secondary particles entering the water, allowing 
the shower direction and thus the direction of the primary particle to be reconstructed. The 
detector and the detector simulation used to study its sensitivity are described in detail 
elsewhere ( [Atkins et al. 1999[) . 



The water Cherenkov technique uses water both as the detection medium and to 
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transfer the energy of air shower photons to charged particles via pair production or 
Compton scattering. Consequently, a large fraction of shower particles can be detected, 
leading to a high sensitivity even for showers with primary energies below 1 TeV. 

Milagrito operated with a minimum requirement of 100 hit tubes per event, where the 
discriminator threshold was set so that a photomultiplier signal of ~0.25 photoelectrons 
would fire the discriminator. The direction of the shower plane is determined with 
an iterative least squares (x^) fitter using the measured times and positions of the 
photomultipher tubes. Only tubes with pulses larger than 2 photoelectrons are used in the 
fit, and in subsequent iterations, tubes with large contributions to are removed. The 
resulting angular resolution is a strong function of the number of tubes remaining in the 
final iteration of the fit, Ufu. If no restriction is made on Ufn, Monte Carlo simulations 
indicate that the median space angle between the fitted and the true shower direction is 
about 1.1° for a source at the declination of MrkSOl. 

The optimal cut on Ufu and the optimal bin size for a point source search depend upon 
the observed nju distribution and the angular resolution as a function of nju. Since the 
point spread function of Milagrito is not well characterized by a two-dimensional Gaussian, 
the standard formulae are inappropriate. To estimate the angular resolution as a function 
of Ufit, the detector is divided into two independent, interleaved portions (similar to a 
checkerboard). For each band of njn, the distribution of space angle differences between 
the two portions of the detector are stored. In the absence of systematic effects, these 
distributions can be interpreted as twice the point spread function of the detector for the 
given band of nfu ( |Alexandreas et al. 1992| ). Under the assumption that the point spread 



function for 7-ray showers is identical to that of hadron- induced air showers, one can use 
the above distributions to determine the optimal cut on njn and the optimal size of the 
angular bin. Figure |I] shows the expected significance of a source as a function of angular 
bin size for three different cuts on Ufu. The analysis indicates that requiring nfu > 40 with 
a bin size of radius 1.0°, which on average contains 57% of the source events, is optimal for 
a binned analysis. As shown in Figurej^, for a rather wide range of cuts, the significance of 
an excess depends only weakly on the chosen source bin size. 

As the detector is much smaller than the typical lateral size of a shower, the shower 
core, i.e. the point where the primary particle would have struck the detector had there 
been no atmosphere, is outside the sensitive detector area for a large fraction of showers 
fulfilling the trigger condition. Assuming a differential flux following for the proton 

background and for a typical 7-source, 16 % of the proton showers and 21 % of the 

7-showers triggering the Milagrito detector have their cores on the pond. This leads to a 
broad distribution of detected events with no well defined threshold energy. Monte Carlo 
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simulations using the MrkSOl spectrum given in Equation|l| predict a distribution starting 
at energies as low as 100 GeV, with 90 % of the detected events having an energy in excess 
of 0.8 TeV. The median energy of detected showers depends on the declination 6 and the 
spectral index of the source, and typical values are 3 TeV for 6 = 39.8° (MrkSOl) and 7 TeV 
for 6 = 22.0° (Crab nebula, assuming an spectrum). 

Detector performance is best evaluated by observations of well-known sources. The 
standard candle of VHE astronomy is the Crab nebula. Simulations indicate that the 
expected statistical significance of the excess above background from the Crab nebula 
in Milagrito is too small to be used for testing Milagrito's performance, and indeed no 
significant excess from this source was observed. However, the large average flux of Mrk 501 
during its flaring state in 1997 results in an expected event rate from MrkSOl 3.6 times 
the Crab rate for Milagrito. MrkSOl can therefore be used to measure the sensitivity of 
Milagrito and to test the reliability of the detector simulation. 

3. Results 

Milagrito took data on MrkSOl from February 8, 1997 to May 7, 1998. The effective 
exposure time was about 370 days, with most of the downtime being due to power outages, 
detector maintenance, and upgrades. Milagrito started operation with about 0.9 m of water 
above the tubes. The water level was increased starting in November 1997 to study how 
the sensitivity changes with water depth. The trigger rate was about 300 Hz with 0.9 m of 
water and increased to 340 Hz (400 Hz) at a depth of l.Sm (2.0m). 

The measured rate of 2420 ± 80 reconstructed events per day for 0.9 m water depth in 
a typical bin with 1.0° radius at the same declination as MrkSOl is in good agreement with 
the predicted rate of 2460l9o° events per day from protons. Helium, and CNO nuclei (the 
error accounts for the uncertainty in the measured flux). The contributions from He and 
CNO to this predicted rate are 27 % and 4 %, respectively. 

The isotropic cosmic ray background flux exceeds the 7-signal from MrkSOl by several 
orders of magnitude. The expected background flux in the source bin must be subtracted 
from the measured one in order to obtain the number of excess events from the source. Since 
the background in the source bin depends on its exposure and the detector efficiency in 
local angular coordinates, the background is calculated directly from the data ([Alexandreas 
et al. 1993|) . For each detected event, "fake" events are generated by keeping the local 
zenith and azimuth angles {6, (p) fixed and calculating new values for right ascension using 
the times of 30 events randomly selected from a buffer that spans about 2 hours of data 
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taking. The background level is then calculated from the number of fake events falling into 
the source bin. By using at least 10 fake events per real event, the statistical error on the 
background can be kept sufficiently small. 

Figure]^ shows the significance of the observed signal as a function of right ascension 
and declination in a 6° x 6° region with the MrkSOl position in the center. For each bin, 
the significance is calculated for the area of the circle with radius 1.0° and the bin center as 
the central point, hence neighboring bins are highly correlated. 

At the source position, 918 954 events are observed with an average expected 
background of 915 330 ± 250 events. The excess of 3624 ± 990 events corresponds to a 
significance of 3.7 sigma. We interpret this result as a reconfirmation of Mrk501 as a TeV 
7-ray source during this period. The corresponding excess rate averaged over the lifetime 
of Milagrito is (9.8 ± 2.7) day"\ The excess rate measured between February and October 
1997 can be directly compared to the 7-rate expected using the average fiux measured by 
atmospheric Cherenkov telescopes during this period. Using the fiux given in Equation |l], 
Monte Carlo simulations of a full source transit predict a 7-rate of (12.5 ± 3.8) day^^, which 
is in good agreement with the measured rate during this period of (13.1 ± 4.0) day^^. 

Figure]^ shows excess divided by background for the lifetime of Milagrito. At 
Milagrito's level of sensitivity, the fiux is consistent with being constant in time. 

The analysis was extended to 10 other nearby blazars {z < 0.06) in Milagrito's field of 
view, including Mrk421, but Mrk501 remains the only analyzed source with a significance 
in excess of 3 a. Results from this blazar sample are reported elsewhere ( |Westerhofi" et al" 
[TO]) . 

4. Conclusions and Outlook 

Milagrito, the first TeV air shower detector based on the water Cherenkov technique, 
observed an excess with a statistical significance of 3.7 a from the direction of Mrk501 
between February 1997 and May 1998. The excess is in agreement with expectations based 
on simulations and indicates that the technique is working as anticipated. 

Milagrito served as a prototype for the full Milagro detector. In its final stage, Milagro 
has a size of 60 m x 80 m x 8 m and two layers of photomultiplier tubes, an upper layer with 
450 tubes at a depth of 1.5 m, and an additional layer with 273 tubes at a depth of 6.2 m. 
With its larger effective area and the ability to reject some of the cosmic ray background, 
Milagro will be at least 5 times as sensitive as Milagrito. Data taking began in early 1999. 



- 7- 



This research was supported in part by the National Science Foundation, the U.S. 

Department of Energy Office of High Energy Physics, the U.S. Department of Energy 
Office of Nuclear Physics, Los Alamos National Laboratory, the University of California, 
the Institute of Geophysics and Planetary Physics, The Research Corporation, and the 
California Space Institute. 

REFERENCES 

Aharonian, F., et al. 1999, A&A, 342, 69 

Alexandreas, D.E., et al. 1992, Nucl. Instr. Meth. A, 311, 350 

Alexandreas, D.E., et al. 1993, Nucl. Instr. Meth. A, 328, 570 

Atkins, R., et al. 1999, to be submitted to Nucl. Instr. Meth. A 

Catanese, M., et al. 1998, ApJ, 501, 616 

Chadwick, P.M., et al. 1999, ApJ, 513, 161 

Gaidos, J. A., et al. 1996, Nature, 383, 319 

Hoffman, CM., Sinnis, C, Fleury, P., and Punch, M. 1999, to appear in Reviews of Modern 
Physics, 71 

McCuUough, J.F., et al. 1999, Proc. 26th Int. Cosmic Ray Conf. (Salt Lake City), 2, 369 
Ong, R.A. 1998, Phys. Rep., 305, 93 

Protheroe, R.J., et al. 1997, Proc. 25th Int. Cosmic Ray Conf., Durban, Highhght Session 
on Mrk501, 8, 317 

Punch, M., et al. 1992, Nature, 358, 477 

Quinn, J., et al. 1996, ApJ, 456, L83 

Quinn, J., et al. 1999, ApJ, 518, 693 

Samuelson, F.W., et al. 1998, ApJ, 501, L17 

Westerhoff, S., et al. 1998, Proc. 19th Texas Symposium on Relativistic Astrophysics, Paris 
(France) 



This preprint was prepared with the AAS IM^jX macros v4.0. 



-8- 



o 
c 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



bin radius [°] 



Fig. 1. — The relative significance of a signal from a point source at the position of MrkSOl 
as a function of the source bin radius for various cuts in the number of tubes in the shower 
plane fit. 
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Fig. 2. — The significance for an event excess as a function of right ascension (RA) and 
dechnation (DEC) in a 6°x6° region with the MrkSOl position (RA=253.468°, DEC=39.760° 
(J2000)) in the center. For each bin, the significance is calculated for the area of the circle 
with radius 1.0° and the bin center as the central point, thus neighboring bins are highly 
correlated. The circle indicates the source bin. 
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Fig. 3. — The fractional event excess from MrkSOl as a function of time. The shaded 
area indicates the time period for which no data from atmospheric Cherenkov telescopes is 
available. 



